Quantum kinetic theory for the collapse of Bose-Einstein condensate 
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We show how to apply the quantum kinetic theory in an inhomogeneous system like harmonically 
trapped Bose condensate. We calculate the temperature dependence of the critical number of 
particles by a quantum kinetic theory within the Hartree-Fock approximation and find that there is 
a dramatic increase in the critical number of particles as the condensation point is approached. Our 
results support the experimental result which was obtained well below the condensation temperature. 
The nature of the temperature dependence of the critical number obtained from the quantum kinetic 
theory is the same as obtained from the scaling theory. 
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For the ultracold Bose gas, the interaction is charac- 
terized by s-wave scattering length a s Atomic in- 
teraction as well as the value of scattering length(a s ) can 
be controlled in Feshbach resonance Stability and col- 
lapse of Bose gas with negative scattering length has been 
observed in the clouds of ultracold 7 Li[ij and 85 Rb@,|(|. 
If the interaction is attractive (a s < 0), the gas tends 
to increase the density of the central region of the trap. 
This tendency is opposed by the zero-point energy and 
thermal energy of the atoms. If the number of atoms is 
greater than a critical number(iV c ), the central density 
increases strongly and the zero-point and thermal energy 
are no longer able to avoid the collapse of the gas. 

To calculate the critical number(7V c ) for the collapse of 
Bose-Einstein condensate (BEC) we do not start from the 
free energy but adopt a mere kinetic theory like approach 
based on the energy and the pressure. In the quantum 
kinetic theory we will start from the Hartree-Fock(H-F) 
energy of the system. We will pick the kinetic energy 
term and two body interaction term of the H-F energy. 
The kinetic energy of the particles causes a outward pres- 
sure and the attractive interaction causes a inward pres- 
sure. At the critical number of particles the two pressure 
would be the same. Beyond the critical number of parti- 
cles the inward pressure would be larger than the outward 
pressure and the system would collapse. In this letter we 
shall explicitly calculate N c by a quantum kinetic theory 
within the H-F approximation. We shall explicitly show 
the temperature dependence of N c . 

Before going into the details of the quantum ki- 
netic theory we estimate the critical number by a semi- 
qualitative theory. We consider a system of Bose gas of 
TV number of particles where all the particles are harmon- 
ically trapped with angular frequency lj. Let the mass 
of each particles be m. As the temperature(T) goes to 
zero, all the particles come down to the ground state and 



the system is well described by the ground state wave 
function v E'o( r ) = 



A e 2? in the position(r) space, 



where I = y'S/mw is the length scale of the oscillators. 
For T = 0, the density of the condensed particles is de- 
scribed as 



n (r) =| * (r) 



N 



(1) 



In absence of collision the number density of the ex- 
cited particles isQ 



where At — 



n T (r) 



2-rrh 2 



1 



(2) 



— T , and gs (x) = x+x 2 /2 3 / 2 +x 3 /3 3/2 + 
... is the Bose-Einstein function of a real variable x. 

Let us consider the interaction potential as V^ n t( r ) = 
gS 3 (r), where g = — 47vh a is the coupling constant and 
a = —a s is the absolute value of the s-wave scattering 
length. For a dilute gas we must have | < 1. The typi- 
cal two body interaction energy for N number of particles 
is ~ N 2 g/2l 3 . For this interaction, the gas tends to in- 
crease the density of the central region of the trap. Well 
below the BEC temperature (T c ), i.e. for T — > 0, this 
tendency is resisted by the zero-point energy (~ NHu>) 
of the atoms. In the critical situation the typical oscilla- 
tor energy must be comparable to the typical interaction 
energy. So, at T = 0, we must have N c frui ~ N 2 g/2l 3 . 
From this relation we can write ~ 1 . 

For < T < T c , the typical total energy of the system 
isQ N 4 / 3 huj. So at these temperatures we must have 



Nc /3 hoj - N 2 g/2l 3 . 



From this relation we have 



[i]l/2 > 1. 

However, near T c the length scale of the system isQ 
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h&!k ^ /TV 1 / 6 . So, near the condensation tem- 

perature we must have Nt^hw ~ N 2 g/2L 3 r . From this 
relation we have ^f^ ~ ] 5 > 1. 



2 



In the above semi qualitative argument we got oc 1 
for T = 0, oc [i]Va for < T < T c and £*a oc 
for T = T c . The proportionality constants were deter- 
mined by the scaling theory 0, The scaling theory 
gave the estimated exponents as well as the proportion- 
ality constants as@, ©] 0.671 for T = 0, @ 1.210 ? 



for < T < T c and [8j 2.253 for T = T c . Since the 
exponents were evaluated by a semi qualitative calcula- 
tion, the evaluation of the exponents is independent of 
the scaling theory. It only determines the proportional- 
ity constants. In the proportionality constants there is a 
temperature dependence factor( p-j^sp ) f° r < T < T c 
Q . Now we are interested to know weather this factor is 
independent of the scaling theory. To do so we will eval- 
uate the proportionality constants by a quantum kinetic 
theory. We will see that the quantum kinetic theory will 
also give the correct exponents and the same tempera- 
ture dependent factor of the proportionality constants. 
So, the scaling theory only gives the pre-factors of the 
proportionality constants. The quantum kinetic theory 
will give some different pre-factors of the proportionality 
constants. 

Within the H-F approximation we have the expression 
of energy functional asQ 



E 



+V(r)n (r) + V(r)n T (r) + |n§(r) 
+2gn (r)n T (r) + gn%{r)] 



(3) 



To evaluate the above energy functional of eqn.(3) 
we have to know n (r) and nxlj) for the interacting 
case. But, for the simplicity of the calculation we keep 
the same form of no(r) and riT(r) as in eqn.(l) and 
eqn.(2). The same form of nj<(v) is obtained from the 
semiclassical form of B-E statistics such thatQ fi(p, r) — 
— 5 where p is the momentum of a single 

particle. With this form of statistics we have the total 
number of excited particles as N = ( fc f J c ) 3 C(3). For 
T < T c , the total number of particles in the ground state 
is no — N — Nt = N(l — rf-) 3 - From the variation in 
statistics we evaluate the energy functional of eqn.(3) as 



E(t) = [(l + l)a(t)-c 2 ]N 



where t = T/T c , c x {t) = |(1 - t 3 ) 

1 N 
/2tt ' 



ml 2 

, 3 Af 1/3 C(4)t 4 



(4) 



f 



5" 



t 9/2 



2 [C(3)]V3 

such that S' 



and 



lC(3)l a I 
SiJ = l (ij)3/2(i+j)3/2 ~ 0.6534. 

From the above eqn.(4) we get the kinetic energy of 



z 



0.2 




T(Tc) 



FIG. 1: Plot of the critical number N c in units of - with 
temperature(T) in units of (T c ). Here | = 0.0066 6]. This 
plot follows form eqn.(9). 



the system as 



E k (t) = a(t)N— 
ml-* 



and get the interaction energy of the system as 

h 2 h 2 

E int (t) = c 2 (t)N— ^ = -c 3 {t)N 2 a— 3 
ml ™ 10 



ml c 



(5) 



(6) 



where c 3 (t) = [1 - t 3 ] 2 + ^^M^N-^H 3 / 2 [1 - t 3 ] + 

For < T < T c , the effective volume of the system is 
V = 4vr/ 3 /3. The outward pressure is P out = ~^r- The 
magnitude of the inward pressure is Pin — ^ qv ± ' the 
critical number of particles Pi n would be equal to P ou t- 
From this equality relation we can write the expression 
of the critical number as 

iVc ~ 3ac 3 (t)M^ j [V 
Eqn.(7) gives the expression of the critical number as 



N c a 
I 



2ci(0) 
3c 3 (0) 



1.253 for T = 0. 



(8) 



Now we see that our scaling result with T > is con- 
sistent with the scaling result of [1, Q for T = 0. The 
numerical result for T = is[l(| ^f 2 - = 0.575 and the 
experimental result for T = is ^f- = 0.459. 

To achieve the Bosc-Einstein statistics the necessary 
condition is such thatQ a/ 1 < 1. For < T < T c , con- 
sidering a/ 1 <C 1, we can approximately write c\{t) = 



3 



| [C(3)] 4 / 3 " Similarly, for < T < T c , we can approx- 
imately write c 3 (t) = ^L=[l - < 3 ] 2 . For < T < T c , 
the expression of iV c obtained from eqn.(7) is ^f 2 - = 

2x i 



1C(3)] 4 / 3 



; [l-t 3 ] 2 



. From this relation we can write 



t b 



N r a . I , -i /o 

-2- = 0.779 - 7 7 ■ . 
Z a (1-i 3 ) 3 



for 0<T<T c . (9) 



The above eqn.(9) is represented in FIG. 1. Now we see 
that the nature of the temperature dependence of the 
critical number is the same as obtained from the scaling 
theory [8(. 

For T = T c , the length scale of the system isQ Lr a = 
IN 1 ' 6 . For this case the effective volume of the system is 
V = 4?d 3 iV 1 / 2 /3. For this reason the critical condition 
would be the same as eqn.(7) except the factor V which 



is to be replaced by 



So, for T = T c the expression 



of N r would be obtained from 



3ac 3 (l) AttNc /2 



1/3 



(10) 



Now we see that the critical number increases dramati- 
cally as the condensation point is approached. 



Initially we qualitatively explained the physics of the 
collapse of the attractive atomic Bose gas. Then we semi 
qualitatively estimated N c a/l for the various range of 
temperatures (0 < T < T c ). Finally we calculate N c a/l 
by a quantum kinetic theory within the Hartree-Fock ap- 
proximation. Our calculation supports the semi qualita- 
tive estimation of N c a/l. The quantum kinetic theory 
gives the same nature of the temperature dependence 
of the critical number as obtained from the scaling the- 
ory. Since - is a very large number and the nature of 
the temperature dependence of the critical number is the 
same for the two different theories, the two theories are 
equal a priory and self consistent with respect to the ba- 
sic physics of the collapse of the Bose gas. However, the 
calculations of the scaling theory is more rigorous than 
that of the quantum kinetic theory. The scaling theory 
result at T = is also closer to the experimental result. 



Putting the value of ci(l) and 03(1) into the eqn.(10) we 
get 



^ = 96.062[-] 5 for T = T C 
I a 



(11) 
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